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ABSTRACT 

The formation of C=C and C=C bonds from the processing of pure c-C(,Hn (cyclo- 
hexane) and mixed f^OrNH^x-CeH^ (1:0.3:0.7) ices by highly-charged, and energetic ions 
(219 MeV 16 7+ and 632 MeV 58 Ni 24+ ) is studied. The experiments simulate the physical 
chemistry induced by medium-mass and heavy-ion cosmic rays in interstellar ices analogs. 
The measurements were performed inside a high vacuum chamber at the heavy-ion acceler- 
ator GANIL (Grand Accelerateur National dTons Lourds) in Caen, France. The gas samples 
were deposited onto a polished Csl substrate previously cooled to 13 K. In-situ analysis was 
performed by a Fourier transform infrared (FTIR) spectrometry at different ion fiuences. Dis- 
sociation cross section of cyclohexane and its half-life in astrophysical environments were 
determined. A comparison between spectra of bombarded ices and young stellar sources in- 
dicates that the initial composition of grains in theses environments should contain a mixture 
of H2O, NH3, CO (or CO2), simple alkanes, and CH3OH. Several species containing double 
or triple bounds were identified in the radiochemical products, such as hexene, cyclohexene, 
benzene, OCN , CO, CO2, as well as several aliphatic and aromatic alkenes and alkynes. 
The results suggest an alternative scenario for the production of unsaturated hydrocarbons 
and possibly aromatic rings (via dehydrogenation processes) in interstellar ices induced by 
cosmic ray bombardment. 

Key words: methods: laboratory - ISM: cosmic rays - ISM: molecules - molecular data - 
astrochemistry - astrobiology 



1 INTRODUCTION 

A large number of gas-phase molecules containing up to 13 atoms 
has been discovered in interstellar space (see review in Herbst & 
van Dishoeck 2009; Ehrenfreund & Charnley 2000). These inter- 
stellar molecules are mainly organic compounds involving unsatu- 
rated carbon-chain species such as simple alkenes (C„H 2 „), cumu- 
lene carbenes (H2C(=C)„), and alkynes (C„H2„-2) (e.g. Thaddeus 
& McCarthy 201 1; Herbst 1995; Winnewisser & Herbst 1987). Cu- 
mulene carbenes and other unsaturated carbon-chain species have 
also been suggested in literature as possible carriers for diffuse in- 
terstellar bands (DIBs) (Apponi et al. 2000). 

Unsaturated small compounds such C4H2, Cgft, 
methylpolyynes (e.g. CH3C2H and CH3C4H) and also ring- 
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compounds such as benzene have also been detected in space 
environments such as in dust shells of late stars, proto-planetary 
nebulaes (Cernicharo et al. 2001a, 2001b). Long aliphatic unsat- 
urated chains containing nitrogen atoms, such as polyacetylenes 
or polyynes H-(C=C)„-H, and cyanopolyynes (H-(C=C)„-CN), 
may also be present in interstellar medium (Fukuzawa et al. 1998). 
Polyacetylenes are considered the missing link between small gas 
phase molecules, among them acetylene (C2H2), carbon monoxide 
(CO), and hydrogen cyanide (HCN), and carbonaceous grain 
particles formed in the outflow of carbon-rich stars (Cordiner & 
Millar 2009; Duley et al. 2005). Aromatic compounds such as 
naphthalene and anthracene cation have been detected in cold 
molecular clouds (Iglesias-Groth et al. 2008, 2010). In addition, 
Duley et al. (2005) suggested the presence of cyclic alkanes and 
alkenes, such as cyclohexene, in attempt to justify the appearance 
of spectral features near 3.4 fim in the infrared spectra of inter- 
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stellar clouds (IC) and in diffuse interstellar medium (DISM). For 
example, they found that the average CH 2 /CH 3 ratio in a number 
of cyclohexene and cyclohexadiene compounds with spectra that 
match the 3.4 yum band in the proto-planetary nebula (PPN) CRL 
618 is «1,15. Single ring systems, such as C-C3H and C-C3H2, have 
also been detected in space (Smith 1992). 

Several mechanisms have been proposed to explain the pres- 
ence of these unsaturated compounds in space including gas-phase 
ion-molecule reactions and neutral-neutral reactions (e.g. Cherch- 
neff et al. 1992; Taylor & Duley 1997; Fukuzawa et al. 1998; Du- 
ley et al. 2005, and references therein). The present experimental 
radiochemical study suggests an alternative scenario for the pro- 
duction of unsaturated carbon chain species (and dehydrogena- 
tion) in interstellar ices induced by cosmic rays bombardment. The 
experiments were performed bombarding cyclohexane (c-C 6 H 12 )- 
containing ices at 13 K with fast ions to simulate the physical chem- 
istry induced by cosmic rays in interstellar unsaturated ice analogs 
inside dense regions of interstellar medium, such as molecular 
clouds and protostellar disks. Although no large saturated cyclic 
compound in the interstellar medium, such as cyclohexane, has 
been direct detected, the presence of this compound, as well as of 
other large saturated hydrocarbons has been suggested in these re- 
gions (e.g. Pendleton et al. 1994; Duley et al. 2005). In this paper, 
this compound is just a prototype for saturated hydrocarbons. We 
expect that the unsaturation processes described here also occur 
with other saturated hydrocarbons already detected in space(e.g. 
Pendleton & Allamandola 2002; Duley et al. 2006). In addition, sat- 
urated hydrocarbons can be formed on the surface of frozen grain 
by H-atom reaction on graphite grains (Bar-Num et al. 1980) and 
direct hydrogenation (catalytic hydrogenation) of alkenes/alkynes 
(e.g. Olah & Molnar 2003; Marcelino et al. 2007). 

Section 2 describes briefly the experimental setup. The results 
on the radiolysis of saturated hydrocarbon containing ices, includ- 
ing formation and dissociation cross sections, are presented and 
discussed in section 3. Some astrophysical implications, as well as 
an estimate for half-lives of these molecules at possible astrophys- 
ical environments, are provided and discussed in section 4. Finally, 
section 5 contains the final remarks and conclusions. 



2 EXPERIMENTAL 

To simulate the physico-chemical changes induced by medium- 
mass and heavy-ion cosmic rays in interstellar ices analogs, we 
used the facilities of the heavy-ion accelerator GANIL (Grand Ac- 
celerateur National dTons Lourds) in Caen, France. The measure- 
ments were performed inside a high vacuum chamber which could 
be coupled to different ion beamlines: IRRSUD (for projectiles 
with 219 MeV 16 7+ ) and SME (for 632 MeV 58 Ni 24+ ). The beam 
flux was (j> = 2 x 10' ions cm" 2 s -1 . The gas samples were de- 
posited onto a polished Csl substrate previously cooled to 13 K. 
The ion projectiles impinged perpendicularly onto the ice target. In- 
situ analysis was performed by a Fourier transform infrared (FTIR) 
spectrometry at different ion fluences. 

In this paper, results for pure C-C6H12 ice and mixed 
H 2 0:NH3:c-C6H 12 (1:0.3:0.7) ice, both at 13 K are presented. Pure 
c-C 6 H l2 ice was bombarded by 219 MeV 16 7+ (13.7 MeV/u). For 
the mixed ice, the projectile employed was 632 MeV 58 Ni 24+ (10.9 
MeV/u). The equilibrium charge state of O and Ni atoms (indepen- 
dent of the initial charge state) after several collisions with matter 
(ice) is around 7 and 13, respectively (e.g. Nastasi et al. 1996). 

In-situ infrared spectra were recorded for ices irradiated at dif- 
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Figure 1. Schematic diagram of the experimental set-up. The ion beam im- 
pinges perpendicularly on the thin ice film deposited on a Csl crystal (taken 
from Pilling et al. 2010a). 



Table 1. Infrared absorption coefficients (band strengths) used in the col- 
umn density calculations for the observed molecules. 



Wavenumber Wavelength 


Assignment 


Band strength 


Ref. 


(cm" 1 ) 


(pm) 




(cm molec ) 




2340 


4.27 


C0 2 (v 3 ) 


7.6xl0~ 17 


[1] 


-2165 


4.62 


OCN~ (v 3 ) 


~4xl0- 17 


m 


2135 


4.67 


CO (v,) 


l.lxl0~ 17 


[31 


1299 


7.69 


CH4 (v 4 ) 


7xKr 18 


[4] 


1095 


9.13 


NH3 (v 2 ; umbrella mode) 


1.2xl0-' 7 


[5] 


861 


11.6 


c-C 6 Hi 2 (-C-C- stretch) 


8X10- 19 


[6] 


-800 


12.5 


H 2 (yi\ libration mode) 


2.6xl0~ 17 


m 


719 


13.9 


AlkeneS (-C-H out-of-plane bending) 


~8xl0~ 18 


[71 



[1] Gerakines et al. 1995; [2] Average value from d'Hendecourt et al. 
(1986) and Demyk et al. (1998); [3] Jiang et al. 1975; [4] Gerakines et al. 
2005; [5] Kerkhof et al. 1999; [6] d'Hendecourt & Allamandola (1986); 
[7] Tentative average value for simple alkenes, such as cyclohexene (see 
details in text). 



ferent fluences, up to 6 x 10 13 ions ctrr 2 using a Nicolet Fourier- 
transformed infrared spectrometer (Magna 550) from 4000 to 600 
cm" 1 with 1 cm" 1 resolution. A background allowing the ab- 
sorbance measurements was collected before gas deposition. Ex- 
perimental details are given elsewhere (Seperuelo Duarte et al. 
2009; 2010; Pilling et al. 2010a; 2010b). 

The sample-cry ostat system can be rotated over 180° and fixed 
at three different positions to allow: i) gas deposition, ii) FTIR 
measurement, and iii) perpendicular irradiation as shown in Fig.Q] 
(taken from Pilling et al. 2010a). The thin ice film was prepared 
by condensation of gases (purity superior to 99%) onto a Csl sub- 
strate attached to a closed-cycle helium cryostat, cooled to 12-13 
K. Experiments without sample bombardment did not show any 
chemical alteration of ice sample due to vacuum with exception 
of a small water layering as a function of time. During the exper- 
iment the chamber pressure was around 2 x 10~ 8 mbar. After the 
irradiation the samples were heated (~ 2 K min~') to room tem- 
perature and some IR spectra was taken during this process. Future 
experiments will be performed to increase the amount of produced 
organic residues which will be analyzed further by ex-situ chro- 
matographic techniques. 

The molecular column densities of samples were determined 
from the relation between optical depth t v = ln(/ //) and band 
strength, A (cm molec -1 ), of the respective sample vibrational 
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Figure 2. FTIR spectra from 4000-600 crrT 1 of non-irradiated ices at 13 K: 
mixed H 2 0:NH 3 :c-C 6 Hi 2 (1:0.3:0.7) ice (top), pure c-C 6 Hi 2 ice (bottom). 
For comparison, spectra of pure H2O ice and pure NH 3 ice are also shown. 

mode. In this expression, I and / are the intensity of light at a 
specific wavenumber, before and after passing through a sample, 
respectively. Because the absorbance measured by the FTIR spec- 
trometer is A v = log(/ //), the molecular column density of ice 
samples is given by 

1 r 2.3 f 

N = — I Tydv = — I a v dv [molec cm -2 ], (1) 
A J A J 

where a v = ln(/ //)/ln(10) = t v /2.3. 

From these measurements and assuming an average density 
for the ice samples of about 1 g/cm 3 , the thickness, the deposition 
rate and the sample mass were determined (see equations in Pilling 
et al. 201 1). For the pure c-CgHii, the thickness was about 1.8 /jm, 
the deposition rate about 1 1 yum/h, and the mass about 100 yug. For 
the mixed H 2 0:NH 3 :c-C6H 12 (1:0.3:0.7)), the estimated thickness 
was 4.5 //m, the deposition rate ~ 27 yum/h, and the mass ~ 250 yug. 

Knowing the ice thickness and the initial column density and 
considering that each monolayer has roughly N M l ~ 10 15 molec 
cnT 2 , the thickness of a single monolayer, d M t, was estimated by 
the expression 

id | d 

d ML = 10 4 = 10 19 — [A], (2) 

NJN ML N„ 

where d is the ice thickness (pan) and N is the initial molecular 
column density (molec cm -2 ). The monolayer thickness of pure 
cyclohexane ice and mixed H 2 0:NH 3 :c-C 6 H 12 (1:0.3:0.7) ice were 
around 15 and 8 A, respectively. 

The analyzed ice layers were thin enough: i) to avoid satura- 
tion of the FTIR signal in transmission mode, and ii) to be fully 
crossed by an ion beam with velocity approximately constant. This 
latter point is important because the relatively low total kinetic en- 
ergy loss of the projectile in the film guarantees that the studied 
cross sections remain constant. 

The dissociation cross sections determined from such exper- 
iments employing ions have a good reproducibility as we observe 
for the two H 2 0:C0 2 (1 : 1) 13 K ices irradiated by 52 MeV Ni ions 
with the same instrumentation (Pilling et al. 2011). For pure ices, 
the temperature is the main parameter that can affect experimental 
reproducibility. However, for mixed ices, besides the temperature, 
changes in the initial molecular abundance ratio may affect the en- 
tire system, resulting in variations in the dissociation cross sections, 



formation cross section of newly produced species, as well as, in 
the sputtering yield (e.g Pilling et al. 2010b). 

For comparison, the infrared (IR) spectra of non-bombarded 
mixed H 2 0:NH 3 :c-C6H 12 (1:0.3:0.7) ice are presented together 
with spectra of non-bombarded pure H 2 0, NH3, and c-CeHi 2 ices 
in Fig. [2] Despite the similarities in the IR spectra corresponding 
to mixed and pure ices, the spectra of mixed ice cannot be ob- 
tained by the simple sum of the spectra of pure ices. The presence 
of different chemical environments in the ice (in the case of mixed 
ices) changes some IR band profiles and also may change the band 
strength of some molecular vibration modes. Except for the sharp 
ammonia peak at 3400 cm -1 (-NH stretching mode), all other sharp 
peaks in the mixed ices are assigned to c-CeH^. The NH 3 um- 
brella vibration mode is easily observed in the mixed ice. How- 
ever, this band is lightly shifted to higher frequencies (1095 cm -1 ) 
in comparison with the band's location observed in pure NH 3 ice 
(1070 cm -1 ). The broad water bands at ~3300 cm -1 (-OH stretch- 
ing mode) and at ~ 800 cm -1 (-OH libration mode) are also easily 
observed in the IR spectra of mixed ices. 

The vibrational band positions and infrared absorption coef- 
ficients (band strengths) used in this work to derive the column 
densities of the selected species are given in Table [T] Because of 
the difficulty in determining the area of NH 3 and H 2 bands due 
to the convolution with other peaks, their column density values 
should be employed with caution. In attempt to quantify the amount 
of produced alkenes, we considered the cleanest and well defined 
alkene peak in the IR spectra, which occurred at 719 cm -1 (13.9 
yum) (=CH out-of-plane bending). For this IR transition, an average 
band strength value of ~8xl0 -18 cm molec -1 was adopted. This 
value is only an approximation for the amount of alkenes, since it 
represents an average value for this band in simple alkenes, such as 
hexene, cyclohexene, benzene (e.g. Pavia et al. 2009; Ruiterkamp 
et al. 2005; d'Hendecourt & Allamandola, 1986). 



3 RESULTS AND DISCUSSION 

Figures|3^-b present the evolution of the IR spectra of pure c-C 6 Hi 2 
ice at 13 K as a function of ion fluence. In both figures, the upper- 
most curve shows the IR spectra of pure c-C6Hi 2 ice at 13 K before 
irradiation. Other curves show the IR spectra at different irradia- 
tion fiuences employing 219 MeV O ions. Each spectrum is offset 
for clearer visualization. The figure inset shows the IR spectrum of 
the pure c-CeHi 2 ice as a whole. The infrared absorption band em- 
ployed to determine the molecular column density of cyclohexane 
is indicated by the arrow. The lowest curve was obtained after a 
fluence of 6 x 10 13 O ions cm . The residual gas inside the cham- 
ber (mainly water) causes the condensation of water (broad bands 
around 3300 cm - ' and 800 cm -1 ), better seen at high fiuences, as 
well as some radiolysis products such as CO (2133 cm -1 ) and C0 2 
(2341 cm -1 ). In both figures, the newly formed products from the 
radiolysis of c-C6Hi 2 are indicated by asterisks and their possible 
assignments are listed in Table lAll 

Figures [4^-b present the infrared spectra of mixed 
H 2 0:NH 3 :c-C 6 H 12 (1:0.3:0.7) ice at 13 K before (uppermost 
curve) and after different irradiation fiuences employing 632 
MeV Ni ions. Each spectrum is offset for clearer visualization. 
As in the previous figure, the infrared absorption band employed 
to determine the molecular column density of cyclohexane is 
indicated by the arrow. The lowest curve was obtained after a 
fluence of 3 x 10 13 Ni ions cm . Figure inset shows IR spectrum 
of non-bombarded c-CgHi 2 ice from 4000 to 600 cm -1 . The broad 
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Figure 3. a) Expanded view of the infrared spectra of pure C-C6H12 ice at 13 
K before (highest curve) and after different irradiation fiuences employing 
219 MeV O ions (medium-mass cosmic ray analog) from 3200 to 1500 
cm -1 . Figure inset shows IR spectrum of non-bombarded C-C6H12 ice from 
4000 to 600 cm -1 , b) Expanded view from 1600 to 600 cm - ' . Ion fiuences 
are indicated. Asterisks indicate the location of newly produced IR bands. 

structure from 3100 to 3500 cm 1 represents a combination of 
vibration modes of water (vi) and ammonia (vi). The narrow peak 
at 2100 cirT 1 is the CO stretching mode (vi). The band at 1600 
cm" 1 consists of two lines corresponding to the water V2 vibration 
mode (1650 cm"') and ammonia v 4 vibration mode (1630 cm" 1 ). 
The feature around 1100 cm" 1 is the umbrella vibration (V2) mode 
of ammonia and the one at 800 cm" 1 is the libration mode (v/J of 
water molecules. The newly formed products from the radiolysis 
of mixed H 2 0:NH 3 :c-C 6 Hi2 (1:0.3:0.7) ice at 13 K are indicated 
by asterisks and their assignments are listed in Table I All 

The formation of C0 2 (2341 cm" 1 ) from the bombardment 
of pure C-C6H12 is sensitively triggered after employing 10 13 ions 
cm" 1 in the ice. In addition, for higher fiuences, this CO2 peak 
has a shoulder in the lower wavenumber region (2334 cm" 1 ) which 
may be attributed to the -C=C- symmetric stretch of non symmetric 
alkynes (Pavia et al. 2009). 

Figuresf5^-b present a comparison of the IR spectra of the two 
irradiated ices (at highest fiuences) in this study with the IR spectra 
of different non irradiated cyclic and aliphatic hydrocarbons from 
NIST databas^] (liquid-phase). Skeletal formulas of each species 
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Figure 4. a) Expanded view of the infrared spectra of mixed H20:NH3:c- 
C6H12 (1:0.3:0.7) ice at 13 K before (highest curve) and after different 
irradiation fiuences employing 632 MeV Ni ions (heavy cosmic ray ana- 
log) from 3200 to 1500 cm" 1 . Figure inset shows IR spectrum of non- 
bombarded mixed ice from 4000 to 600 cm" 1 . b) Expanded view from 1600 
to 600 cm" 1 . Ion fiuences are indicated. Asterisks indicate the location of 
newly produced IR bands. 



(cyclohexene, 1-3 cyclohexadiene, benzene, hexane, 1-hexene, and 
1-hexyne) are also shown. Asterisks indicate the peaks that have 
possible identification in the spectra of irradiated ices. Figure [5^ 
shows the IR spectra from 3200-1500 cm" 1 and Fig.[5p presents the 
1500-600 cm" 1 wavenumber range. The presence of a broad struc- 
ture at 1600 cm" 1 indicates the presence of the rings of aromatic 
systems in the ice after bombardment as result of dehydrogenation 
of cyclohexane molecules. A theoretical study about the energies 
and the transient states involved during the ring opening of cyclo- 
hexane to produce 1-hexene (c-C6H 12 — > 1-C6H 12 ) was performed 
by Sirjean et al. (2006). The authors also suggested that species 
such as birradicals (»C6Hi2» and «C 4 H8») and ethylene (C2H4) are 
present in such ring open mechanisms. This last species is required 
in the production of aromatic compounds in interstellar environ- 
ments (e.g. Frenklach & Feigelson 1989). 

Table I A 1 1 lists the new IR bands formed after ion bom- 
bardment of the cyclohexane containing ices at 13 K. A large 
number of vibrational modes are attributed to alkenes and pos- 
sibly to ring/aromatic systems. For the mixed H20:NH3:c-C6Hi2 
(1:0.3:0.7) ice, the formation of CO, C0 2 , and OCN" are observed. 
Due to the local heat and energy delivered by projectiles in the sam- 
ple, the bombardment is expected to promote desorption of molec- 
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Figure 5. Comparison between the IR spectra of the two irradiated ices at highest fluences (this work) with IR spectra of different non irradiated cyclic and 
aliphatic hydrocarbons from the NIST database. Skeletal formulas of each species (cyclohexene, 1-3 cyclohexadiene, benzene, hexane, 1-hexene, 1-hexyne) 
are shown. Asterisks indicate the peaks that have possible identification in the spectra of irradiated ices, a) From 3200 to 1500 citT 1 . b) From 1600 to 600 



ular hydrogen, however no measurement of this species has been 
performed. 



3.1 Cross section and radiolysis yield 

The column density variation as a function of ion fluence for each 
condensed molecular species (parental or daughter) during radiol- 
ysis can be written as 



dN > V AJ , 



cr dj N, - Y,Q.,{F) 



(3) 



where 2* crfjtNk represents the total molecular production rate of 
the i species directly from the k species, L, is the layering, tTjj the 
dissociation cross section, F, the sputtering yield, and 0,(F) the 
relative area covered by the i species on the ice surface. For more 
details see equations in Pilling et al. (2010a) and de Barros et al. 
(2011). 

For pure cyclohexane ice, the solution of Eq. [3] considering 
no extra production(<x/,i = 0), full covering (flt(F) = 1) and, no 
layering (L=0), is given analytically by 



N = (N + — )exp(-cr,,F)- — 



(4) 



where N , Y, and <x rf are the initial column density, sputtering and 
dissociation cross section of cyclohexane species, respectively. 



In the case of mixed H 2 0:NH 3 :c-C6H 12 (1:0.3:0.7) ice, after 
the rapid ice compaction phase at the beginning of the irradiation, 
water layering tends to recover the ice surface (fl,(F) — > 0) pro- 
gressively preventing the sputtering of other species. Therefore, un- 
der the assumption that both sputtering and layering are negligible 
(Y=L=0) the solution of Eq.[3]for the column density of c-CeHi2 
as a function of fluence is obtained analytically by 



N = N exp(-cr d F) 



(5) 



where N , and <r d are the initial column density, and dissociation 
cross section of cyclohexane species in the mixed ice, respectively. 

For water (and viscous liquids), it may occur an additional 
layering due to the molecular retention in the chamber walls. In 
this case, and considering no other formation source (tr/,a = 0), 
the column density evolution with the fluence is obtained directly 
by solving Eq.[3] which gives 



N = (N - Noo) exp(-cr rf F) + N K 



(6) 



where No, and cr rf are the initial column density and dissociation 
cross section of H 2 0, respectively. In this equation, N m = (L - 
Y)j(Ti is the asymptotic value of column density of water due to 
the presence of layering. In experiments where high layering values 
are observed, the sputtering may be considered negligible. 

As discussed by de Barros et al. (2011), for low fluences (up 
to ~ 5 x 10 12 ions cm -2 ), the column density evolution of daughter 
species produced by the radiolysis of both pure and mixed ices is 
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Table 2. Stopping power and penetration depth values for 219 MeV (13.6 
MeV/u) 16 ions and 632 MeV (10.8 MeV/u) 58 Ni ions in the studied ices, 
calculated by the SRIM code. For comparison, values for 632 MeV (10.8 
MeV/u) 58 Ni, 1 MeV protons, and 12 MeV protons all in pure cyclohexane 
ice, and 1 MeV protons in pure water ice are also shown. 



Projectile 




Stopping 


power 


Penetration 




Electronic 


Nuclear 


Total 


depth (jum) 




(keV/jum) 


(keV/pm) 


(eV/(molec/cm 2 )) 




219 MeV 0" 


250.4 


0.13 


3.49 xlO" 13 


509.7 


632 MeV Ni 6 


2787 


1.70 


1.27 xl0~ 12 


172.5 


632 MeV Ni a 


2802 


1.72 


3.88 xlO" 12 


169.8 


1 MeV p a 


28.89 


0.023 


3.12 xl0~ 15 


21.4 


1 MeV p c 


25.11 


0.021 


7.51 xl0~ 13 


25.1 


12 MeV p" 


4.296 


0.002 


4.03 xl0~ 14 


1540 



Employed ice density = 1 .0 g cm 3 . 

" in pure c-CgHi2 ice (~ 7 xlO 21 molec cm -3 ). 

b in mixed H 2 0:NH 3 :c-C 6 H 12 (1:0.3:0.7) ice (~ 2 XlO 22 molec cm" 3 ). 
c in pure H2O ice (~ 3 XlO 22 molec cm -3 ). 



N k * N o- f , k (F - H±p* F *} (7) 

where Nk is the column density of daughter species k at a given 
fluence. cry * an d (T d,k are the formation cross section and dissocia- 
tion cross section of daughter species k, respectively. No and a d are 
the initial column density and dissociation cross section of parental 
species, respectively. This equation assumes a negligible value for 
the sputtering of daughter species (5^=0). 

Figure [6] shows the evolution of the column density of cyclo- 
hexane in pure and mixed studied ices as a function of ion fluences. 
The column densities of four radiolysis products (CO, CO2, CH4, 
and OCN) observed in the radiolysis of the mixed ice are also 
shown. The decreasing of the cyclohexane column density is re- 
lated to the formation of other species and to the sputtering induced 
by heavy-ions (Seperuelo Duarte et al. 2009; 2010; Pilling et al. 
2010a; 2010b). 

Table [2] shows the stopping power and penetration depth val- 
ues for 219 MeV (13.6 MeV/u) lfi O ions and 632 MeV (10.8 
MeV/u) 58 Ni ions in the studied ices, calculated by the Stopping 
and Range of Ions in Matter - SRIM code0. The SRIM code is a 
collection of software packages that calculate many features of the 
transport of ions in matter (Ziegler et al. 2011). For comparison, 
values for 632 MeV (10.8 MeV/u) 58 Ni, 1 MeV protons and 12 
MeV protons, all in pure cyclohexane ice, and 1 MeV protons in 
pure water ice are also shown. The energy delivered per micron by 
a single 632 MeV Ni ion is roughly 1 1 times higher than for a sin- 
gle 219 MeV O ion, and about 1 10 times higher than for a single 1 
MeV proton. 

The radiochemical formation yield (G/) of a given compound 
per 100 eV of deposited energy, at normal incidence, is written by 

G f = 100y molecules per 100 eV. (8) 

where cry is the formation cross section and S is the stopping power, 
in units of eV/(molec/cm 2 ) (Loeffler et al. (2005). 

This definition can be extended to the radiochemical dissoci- 
ation (destruction) yield of a given compound per 100 eV of de- 
posited energy (Gj) by replacing the formation cross section in 

2 http://www.srim.org/ 



Eq.[8]with the negative value of the dissociation cross section (- 
(Tj). Therefore, negative Gj values indicate that molecules are be- 
ing dissociated or destroyed after energy deposition into the ice. By 
adopting the 5 values from the stopping and ranges module of the 
SRIM code, the values of the radiation yield G in the experiments 
can be determined and compared with values in the literature. 

Figures |6^-b present the column density evolution of species 
initially present in pure and mixed ice experiments, as well some 
newly formed products due to the radiolysis, as a function of ion 
fluences. Fig. [6ji shows the evolution of cyclohexane, the newly 
formed CH4, as well as the possible alkenes produced from the 
bombardment of the pure cyclohexane ice by 219 MeV O ions. The 
evolution of cyclohexane, water and ammonia, four newly formed 
species (CO, CO2, CH4 and OCN~), and possible alkenes observed 
in the radiolysis of the mixed ice by 632 MeV Ni ions are given 
in Fig. |6jj. The lines indicate the fittings using Eq. [4] (pure cyclo- 
hexane), Eq. [5] ( cyclohexane and ammonia in mixed ice), Eq. [6] 
(water in mixed ice) and, Eq.[7](products). Numeric labels indicate 
the models employed in which the parameters are listed in Table[3] 

Table [3] list the values obtained by best-fitting curve employ- 
ing Eqs.|4]|5]|6]and,|7]to the column density evolution data of bom- 
barded ices. First columns show the selected species, followed by 
the formation cross section (only for daughter species), dissocia- 
tion cross section and, radiolysis yields (formation and destruc- 
tion). For the radiolysis yields we employed the stopping power 
values listed in Table|2] 3.49 x 10~ 13 eV/(molec/cm 2 ) (for pure ice) 
and 1.27 x 10~ 12 eV/(molec/cm 2 ) (for mixed ice). In addition, for 
mixed ice it was considered that the projectile energy was absorbed 
by the whole bulk sample and not only for a single species. The last 
three columns in Table[3]are the sputtering yield (only for pure ice), 
the layering and the initial column density. The sputtering value for 
water of 10 4 molec per impact was taken from Brown et al. (1984). 
The water layering was derived by the equation L = N m (Td + Y- 

The destruction cross section obtained in this study (around 
10" 13 cm 2 ) for the frozen species bombarded with 632 MeV Ni ions 
are in the same order of magnitude of the values obtained previ- 
ously also employing swift heavy-ions (Pilling et al. 2010a; 2010b; 
2011). However, the formation cross sections and also the forma- 
tion yields of new species decrease with the projectile energy. This 
issue reinforces that an optimal energy associated with the forma- 
tion of new (organic) species during the processing of astrophysical 
ices could exist (see also Pilling et al. 201 1). 

The formation cross section of alkenes, considering the C- 
H out-of-plane bending mode at 719 cm [ (A ~ 8 x 10~ 18 cm 
molec 1 ), in the pure c-CeH^ ice and mixed H20:NH3X-C 6 Hi2 
ice were about 1 x 10~ 15 and 1 x 10~ 14 cm 2 , respectively. These can 
be compared with similar experiments employing other saturated 
hydrocarbons such as ethane, propane and cyclopropane. Future 
experiments will compare the formation cross section of alkenes 
(-C=C-), induced by cosmic rays in interstellar ice analogs, with 
respect to the number of carbons of parental saturated hydrocarbon 
present initially in the ice. 



4 ASTROPHYSICAL IMPLICATIONS 

Despite the ion flux employed in laboratory experiments be several 
orders of magnitude higher than the flux of similar ions in space 
and, the thickness of laboratory ice be also higher, the energy de- 
livered and damage induced by similar ions in both scenarios are 
similar. For example, lets quantify the role of ice size in the experi- 
ments. As discussed in section 2, the thickness of pure cyclohexane 
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Table 3. Cross sections (ay and crj), radiochemical yield (G/ and G^), and sputtering yield (Y) for cyclohexane (and H2O, NH3 for the mixed ice) and selected 
daughter species obtained from ion bombardment experiments of 13 K ices. The layering (L) and the initial column density (No) are also listed. The fittings 
are shown in Fig. [6] 

Pure c-C 6 Hi2 ice irradiated with 219 MeV 7+ 
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Mixed H20:NH 3 :c-C6H 12 (1:0.3:0.7) ice irradiated with 632 MeV Ni 2 



species" 


°7 

(lO^ 13 cm 2 ) 


(10~ 13 cm 2 ) 


a, 

1 molec 1 
M00 eW> 


-Gd 
1 molec ^ 
M00eV ; 


Y 

qq4 molec n 
v ion ' 


L 

/jq4 molec s 
^ ion ' 


No 

,mis molec ^ 
uu cm 2 > 


Model 


H 2 


0" 


- 3 


0* 


-24 


\ d 


300 


3.3 


4 


NH 3 


0" 


~ 3 


0* 


- 24 








0.86 


5 


c-C 6 H 12 


0" 


- 2 


0* 


- 16 








2.0 


6 


CH 4 


0.14 


0.9 


1.1 


7.1 











7 


CO 


0.07 


< 0.01 


0.55 


<0.08 











8 


OCN" 


0.007 


< 0.01 


0.06 


<0.08 











9 


C0 2 


0.006 


<0.1 


0.05 


<0.8 











10 


Alkenes c 


-0.1 


- 1.5 


-0.8 


- 12 











11 



" Band position given in Table[T] 

* Considering a negligible production of parental species during the radiolysis. 

Considering the C-H out-of-plane bending mode at 719 cirT 1 (A - 8 X 10~ 18 cm molec - '). 
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ice was about 1.5 fim and, the stopping power of the employed pro- 
jectile in this ice was 250 keV /jnT 1 (see Table 2). Thus, the energy 
difference between the projectile at the ice surface (first layer) and 
after goes through the ice was 0,375 MeV (~0,2% of projectile en- 
ergy). For the mixed ice, the thickness was around 4.5 /jm and the 
projectile stopping power was 2788 keV yum -1 , resulting in the pro- 
jectile's energy difference, between top and bottom ice layers, of 
about 12,5 MeV (~2% projectile energy). Therefore, in both cases 
we consider a constant projectile energy within the sample and also 
a linear energy deposit. These assumption are in a good agreement 
with cosmic ray impact in submicron ices in interstellar or inter- 
planetary regions. 

In a similar experiment performed by Pilling et al. (2010a), the 
authors discussed that each projectile induces significant changes in 
the sample only in a region with roughly 3 nm of diameter. Consid- 
ering a constant and homogeneous ion flux of 2x 10 9 ions ctrr 2 s -1 , 
the average distance between two nearby impacts is roughly 300 
nm (about one hundred times higher the length of processed sample 
by a single projectile hit). From these values, we estimate that the 
probability of an ion to hit a given area of 3 nm of diameter in each 
second is about P ~ 1.4 x 10 -4 . Therefore, only after about two 
hours of continuous bombardment we expect a second projectile 
hit in the same nanometric region. This is enough time to consider 
that, in laboratory, each projectile hit always a thermalized region. 
Such scenario is very similar to interstellar or interplanetary low 
flux conditions (over a very extended period of time). 

The equivalent cosmic ray exposure time (Tjsm) in interstellar 
medium with respect to the experimental ion fluence can be ob- 
tained by the equation: 

T ISM - 3.2 x lO" 8 - [year] (9) 



where F is the ion fluence employed in the experiments, in units 
of ions cm 2 , and <f> indicates the flux of galactic cosmic rays (e.g. 
for medium mass cosmic rays 4> M cr or f° r heavy cosmic rays 4>hcr) 
in units of ions cnT 2 s -1 , Pilling et al. (2010a) estimated the inte- 
grated cosmic ray flux for heavy-ion (12 < Z < 29) component with 
energies between 0.1-10 MeV/u in interstellar medium. The value 
found was </>hcr ~ 5 x 10~ 2 cnT 2 s . Pilling et al 201 1, employing 
a similar methodology, estimated the medium-mass (3 < Z < 11) 
component of galactic cosmic rays in interstellar medium with 
the same energy range. The value obtained by those authors was 
<Pmcr ~ 4x 10"' cnT 2 s . Following those authors, both values are 
also a good estimative for the heavy-ion and medium-mass compo- 
nents of galactic cosmic rays at the outer border of the heliopause. 

The chemical evolution of the studied ices extrapolated to 
space conditions (dust grains in interstellar medium) is shown in 
Fig. [7^-b. In this figure, lines are only to guide the eyes. The col- 
umn density ratios of cyclohexane and some of its radiolysis prod- 
ucts over initial cyclohexane column density, as function of equiva- 
lent exposure time in ISM by cosmic rays, are shown in Fig.|7^ for 
pure c-CjH^ ice at 13 K irradiated by 219 MeV O ions (medium- 
mass cosmic ray analog), and in Fig. |7J) for H 2 0:NH3:c-C 6 Hi2 
(1:0.3:0.7) ice at 13 K irradiated by 632 MeV Ni ions (heavy cos- 
mic ray analog). The total column density of produced alkenes was 
obtained by using the area of the IR feature at 719 cm"' (-CH bend- 
ing out plane) considering, as a first hypothesis, that this band can 
be representative of all the alkenes in the ice. This value represents 
an intermediate value between the strongest and weakest IR fea- 
tures in the case of clyclohexane (e.g. d'Hendecourt & Allaman- 
dola, 1986). The error bars in the figure reflect the uncertainty in 
the methodology employed to derive the column density of alkenes 
from IR spectra. The maximum production of alkenes is obtained 
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Figure 6. a) Variation of the column density of cyclohexane and selected 
daughter species (CH4) in pure ice experiment as function of ion (219 MeV 
O) fluence. b) Variation of the column density of cyclohexane, water and 
ammonia and selected daughter species (CO, CO2, CH4 and OCN~) in 
mixed ice experiment as a function of ion (632 MeV Ni) fluence. The lines 
indicate the fittings using Eq.[4](pure cyclohexane), Eq.[5](cyclohexane and 
ammonia in mixed ice), Eq. [6] (water in mixed ice) and, Eq.[7](products). 
The model parameters are given in Table [3] 
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Figure 7. Column density ratios of cyclohexane and some radiolysis prod- 
ucts over initial cyclohexane column density as function of equivalent cos- 
mic ray exposure time at ISM. a) pure c-CgH^ ice at 13 K irradiated by 
219 MeV O ions (medium-mass cosmic ray analog), b) mixed H20:NH 3 :c- 
C 6 Hi2 (1:0.3:0.7) ice at 13 K irradiated by 632 MeV Ni ions (heavy cosmic 
ray analog). In each figure, the half-life of cyclohexane (h/2) in ISM as a 
result of cosmic ray bombardment is indicated. Lines are only to guide the 
eyes. 



after about 3-5 xlO 6 years, independently of the projectile type. 
Its maximum value was roughly the same in both experiments, ~ 
10 2 alkenes per cyclohexane molecule. However, this value must 
be adopted with caution since to the employed assumption for the 
column density of alkenes may not be representative for all alkenes. 
Better determination of the band strengths of alkenes for the 719 
cm -1 IR feature in ice phase is needed to make this methodology 
more accurate. 

Employing the estimated cosmic ray ion flux and the deter- 
mined dissociation cross section, the half-life of cyclohexane (ti/2) 
in ISM as a result of cosmic ray bombardment for pure cCgH^ ice 
is about 5 x 10 6 years (considering only medium-mass cosmic rays) 
and for mixed ice is about 1 x 10 6 years (considering only heavy- 
ion cosmic rays). These values are indicated in Fig. [7] The results 
also show that for the mixed ice (H 2 0:NH 3 :c-C 6 H 12 ), after 20 x 10 6 
years in ISM, almost 20% of the initial cyclohexane was converted 
into CO by heavy cosmic rays, 3% was transformed into OCN~ and 
1% into CO2. This suggests that highly hydrogenated hydrocarbons 
in water-rich grain mantles inside interstellar clouds can be largely 
converted into CO during the lifetime of the cloud. 

Figure |7p also suggests that some daughter species such as 
CH 4 (1300 cm" 1 ) and OCN~ (2165 cirT 1 ) can be used to estimate 
the integrated dose of incoming radiation and, assuming a constant 
cosmic ray flux over the time, the exposure time of interstellar ices 
to the cosmic rays. At ion fluences higher than 3 x 10 12 ions cm 2 



(~ 1 x 10 16 years in ISM), the abundance of these species increases 
almost linearly with the fluence. 



4.1 3.4 fim band 

The 3.4 yum (~2925 cm" 1 ) band observed in the IR spectra of sev- 
eral astrophysical sources, including protostellar ices, hot grains in 
diffuse interstellar medium, dust grains in planetary nebulae and 
meteorites extractions, is frequently associated with hydrocarbons 
(e.g Allamandola et al. 1992, and references therein). As pointed 
out by Pendleton and Allamandola (2002), although the 3.4 yum 
band probes only the CH stretching modes of the carrier, the anal- 
ysis of the peak position, and the band profile has revealed that this 
portion of the refractory material contains aliphatic hydrocarbons 
(open chain -CIT- and -CH 3 groups) that incorporate at least 3% 
of the available interstellar carbon. In addition, the other IR fea- 
tures near this band, for example around 3.2-3.3 fim, also indicate 
the presence of alkenes and aromatic rings (=CH 2 sp 2 ). 

Figure [8]presents a comparison between the 3.4 /jm band ob- 
served in four different astronomical objects: Young Stellar Ob- 
ject (YSO) NGC 7538 IRS9 (Allamandola et al. 1992); Plane- 
tary Nebulae (PN) IRAS 05341+0852 (Joblin et al. 1996); Proto 
Planetary Nebulae (PPN) CRL 618 (Chiar et al. 1998); and Dif- 
fuse Interstellar Medium (DISM) at galactic center through IRS 6E 
(Pendleton et al. 1994), and in the Murchison meteorite extract (de 
Vries et al. 1993) with the FTIR spectra of the studied H 2 0:NH 3 :c- 
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Figure 8. Comparison between the 3.4 /im band of four different astronomical observations (Young Stellar Object NGC 7538 IRS9 (Allamandola et al. 1992); 
Planetary Nebula IRAS 05341+0852 (Joblin et al. 1996); Protoplanetary Nebula CRL 618 (Chiar et al. 1998); Diffuse interstellar medium at galactic center 
IRS 6E (Pendleton et al. 1994), and the Murchison meteorite extract (de Vries et al. 1993) with the FTIR spectra of H 2 0:NH 3 :c-C 6 Hi2 ice (this work) at: a) 
three different ion fluences; b) during sample heating from 13 K up to room temperature. 
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Figure 9. Comparison between the radiolysis products from two simulated 
interstellar ices at 13 K and after warming (H 2 0:NH 3 :CO (1:0.6:0.4) bom- 
barded with 46 MeV 46 Ni 13+ (Pilling et al. 2010a) and H 2 0:NH 3 :c-C6Hi 2 
(1:0.3:0.7) bombarded with 632 MeV 58 Ni 24+ (this work). Each set in this 
figure also compares the laboratories spectra with one astronomical IR spec- 
trum (Young Stellar Object NGC 7538 IRS9 (Allamandola et al. 1992); 
Planetary Nebulae IRAS 05341+0852 (Joblin et al. 1996)); Murchison me- 
teorite extract (de Vries et al. 1993). 

C6H12 ice. Each set in Fig[8ji shows the 3.4 yum band of bombarded 
H 2 0:NH 3 :c-C 6 H 12 ice at three different ion fluences (0, 1.3xl0 13 , 
and 3xl0 13 ions cm -2 ). Each set of Fig.[8p illustrates the evolution 
of 3.4 yum band of bombarded H 2 0:NH 3 :c-C 6 H 12 ice after 3xl0 13 
ions ctrr 2 and at some temperatures during the heating from 13 K 
to room temperature. Temperature values are indicated in the fig- 
ure. 

The 3.46 yum (~2880 ctrr 1 ) sharp feature observed in the YSO 
NGC 7538 IRS9 (squares in the upper set of Fig. [8ji-b) indicates 
that cyclohexane (and possibly other aliphatic hydrocarbons) was 
not present initially in the ices of this protostellar object. As dis- 
cussed by Allamandola et al. (1992), the peak at 3.54 fxm (-2850 



cm -1 ) is attributed to methanol and this species is expected to be in 
the chemical inventory of protostellar grains. 

The comparison between the 3.4 yum band in the PN IRAS 
05341+0852 (the IR spectrum was inverted to appears like an 
absorption-type spectrum in Fig. [8) and the H 2 0:NH3:c-C6Hi 2 ice 
is remarkably good. The -CH 3 stretching mode of aliphatic hydro- 
carbons at 3.37 and 3.48 yum is a good indicator of the amount of 
energy that was delivered to the grains. The same is also true for the 
=CH 2 (sp 2 ) vibration mode at 3.33 /jm (~3010 cm -1 ), which be- 
comes very strong only at higher Fluences (e.g. 3xl0 13 ions cm -2 ). 
As illustrated in Fig.[8p, this band is very sensitive to the tempera- 
ture of the ice, suggesting that the carrier must be a volatile alkene. 
This band may also be used as a thermometer for bombarded inter- 
stellar ices. In addition, since this band is not observed in Murchi- 
son meteorite extraction, this suggests that possible volatile alkene 
species that were present at one time are now completely lost. 
Looking at the temperature effect on the 3.33 /jiti band (4th graph 
set in Fig. [Hp), the interstellar grains in the DISM through IRS 6E 
should have experienced temperatures higher than 50 K. 

In general, the 3.4 yum band in PPN CRL 618, DISM IRS 6E, 
and Murchison extract is very similar with each other, and by ex- 
cepting the sharp absorption peak at 3.37 (~ 2960 cm -1 ), they 
are also very similar to the IR spectra of bombarded H 2 0:NH 3 :c- 
C6Hi 2 ice at a fluence of 3xl0 13 ions cm -2 . The enhancement of 
the 3.37 yt/m band as a function of fluence suggests that these inter- 
stellar grains (and the grains in the Murchison meteorite) probably 
were exposed to a heavy cosmic ray fluence higher than 3xl0 13 
ions cm -2 or to some kind of ionizing agents which delivered an 
energy dose higher than ~ 2 x 10 8 J kg -1 « 10 15 eV ng -1 (10 13 
ions cm -2 ~ 6x 10 7 J kg -1 ; Pilling et al. 201 1). 

Figure [9] presents a comparison of the 3.4 yum band observed 
after heavy-ion bombardment of two simulated interstellar ices 
at 13 K and after warming (H 2 0:NH 3 :CO (1:0.6:0.4) bombarded 
with 46 MeV 46 Ni 13+ (Pilling et al. 2010a) and H 2 0:NH 3 :c-C 6 H 12 
(1:0.3:0.7) bombarded with 632 MeV 58 Ni 24+ (this work). Both ex- 
periments were performed employing the same instrumentation and 
methodology. The ion fluences and ices temperatures are indicated 
in the figure. In each set of this figure we also present a com- 
parison of the laboratories spectra with an astronomical IR spec- 
trum (Young Stellar Object NGC 7538 IRS9 (Allamandola et al. 
1992); Planetary Nebulae IRAS 05341+0852 (Joblin et al. 1996)); 
Murchison meteorite extract (de Vries et al. 1993). The 3.46 yum 
feature observed in the young stellar object NGC 7538 IRS9 is in 
a good agreement with the band detected in the radiolysis prod- 
uct from H 2 0:NH 3 :CO ice at 13 K. However, small features at 3.7 
yum and 3.43 yum seem to match the infrared features of the 300 
K residue for bombarded H 2 0:NH 3 :c-Cf;H 12 ices. Furthermore, an 
additional abundance of methanol is needed to get a good fitting 
of the observation at 3.53 yum. This scenario suggests that the ini- 
tial abundance of this protostelar disk may contain species such as 
NH 3 , CO, and simple aliphatic hydrocarbon including cyclohexane 
and methanol as pointed out by Allamandola et al (1992). The data 
also indicates that warmer ices (closer to the star) may be rich in 
simple aliphatic hydrocarbon and colder ices rich in CO. 

For the planetary nebulae IRAS 05341+0852 (IR spectrum 
was inverted to appears like an absorption-type spectrum in Fig. [9}, 
the best match occurred with the residue (~ 130 K) obtained after 
the ion bombardment of H 2 0:NH 3 :c-C6Hi 2 ice. This puts a con- 



3 Considering a typical density of 1 g cm 3 and a size of 0. 1 yum, the typical 
mass of interstellar grains is about 1 ng. 
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Figure 10. Comparison between IR spectra from 2000 to 1000 cm (6-10 /im) some of interstellar and laboratory ices. The top three curves are infrared 
spectra of young stellar sources obtained by the Infrared Space Observatory (ISO). Lower traces indicate different laboratory spectra of two ammonia- 
containing ices irradiated by heavy-ions at 13 K (H 2 0:NH 3 :CO (1:0.6:0.4) bombarded with 46 MeV 46 Ni 13+ (Pilling et al. 2010a) and H 2 0:NH 3 :c-C6Hi 2 
(1:0.3:0.7) bombarded with 632 MeV 58 Ni 24+ (this work)). 



strain on the temperature of dust grains in the region and indi- 
cates that a negligible amount of CO was in the grains initially. 
The organic residues for both bombarded ices (H 2 0:NH 3 :CO and 
H 2 0:NH 3 :c-C6Hi 2 ) at temperatures hotter than 250 K have some 
similarities in this wavelength region with peak around 3.53 yum 
and 3.42yum associated with aliphatic hydrocarbons. However, the 
IR spectrum of Murchison extraction around 3.4//m can not be fit- 
ted by any spectrum taken from these two simple interstellar ice 
analogs. 



4.2 5-10 /im range 

Figure [TJJ] presents a comparison between IR spectra from 5 to 
10 fjm (2000 to 1000 cnT 1 ) of interstellar and bombarded labora- 
tory ices. The top three curves are infrared spectra of YSOs (W33, 
AFGL989, NGC7538 IRS9) obtained by the Infrared Space Obser- 
vatory (ISO). The four bottom curves in this figure indicate differ- 
ent laboratory spectra of two ammonia-containing ices irradiated 
by heavy and energetic ions: H 2 0:NH 3 :CO (1:0.6:0.4) at fluence 
of 2xl0 13 ion cm 2 (Pilling et al. 2010a), and H 2 0:NH 3 :c-C 6 H 12 
(1:0.3:0.7) at fluence of 3xl0 13 ion cm 2 (this work). For compari- 
son, the organic residues (after heating) of these two simulated ices 
after bombardment are also shown. 

The 7.7 /i m (1300 ctrT 1 ) feature in YSOs, which is attributed 
to CH 4 (-CH deformation mode; Boogert et al. 1997), is best rep- 
resented by the products coming from the radiolysis of ices con- 
taining alkanes initially, rather than by the ices where the primary 
source of carbon comes from CO. This suggests that saturated hy- 
drocarbons, such as cyclohexane, may be a part of the initial chem- 
ical inventory of interstellar ices. However, the presence of broad 
features at 1450 cm"' and 1680 cm" 1 in the astronomical observa- 
tions, similar to the features observed in the IR spectra of the irra- 
diated H 2 0:NH 3 :CO ice (Pilling et al. 2010a) requires the presence 



of CO (or/and C0 2 ) as one of the parental species in the interstellar 
grain ice mantles. Based on these results and with the fact that the 
3.4 fim band in these YSOs requires the presence of methanol (e.g. 
Allamandola et al. 1992), we propose that the initial composition 
of grains in the disc of these young stellar sources should contain a 
mixture of H 2 0, NH 3 , CO or (C0 2 ), alkanes, and CH 3 OH. A future 
work will focus on the initial abundance of these simple compound 
in young stellar environments. 



5 CONCLUSION 

We performed an experimental study of the interaction of medium- 
mass ions (219 MeV 16 O v+ ) and heavy-ions ions (632 MeV 58 Ni 13+ ) 
with cyclohexane-containing ices at 13 K (pure-c-CeH 12 and 
H 2 0:NH 3 :c-C 6 Hi 2 (1:0.3:0.7)) to simulate the physico-chemical 
changes (e.g. molecular unsaturation) induced by cosmic rays in 
saturated hydrocarbon-rich ice analogs inside dense regions of the 
interstellar medium such as dense molecular clouds or protoplane- 
tary disks, as well on the surfaces of solar system ices. Our main 
results and conclusions are the following. 

(i) For the experiment with medium-mass cosmic ray analogs 
(219 MeV O v+ ) on pure c-C 6 H 12 ice at 13 K, the determined dis- 
sociation cross section is 1 xl0~ 14 cm 2 , and the sputtering yield is 
10 3 c-C6H 12 molecules per ion impact. For the experiment in which 
mixed H 2 0:NH 3 :c-C 6 H 12 (1:0.3:0.7) ice at 13 K was bombarded by 
heavy-ion cosmic ray analogs (632 MeV Ni 24+ ), the obtained dis- 
sociation cross sections of H 2 0, NH 3 and c-C 6 H 12 are ~3 xlO -13 , 
~3 xlO" 13 , and ~2 xl0~ 13 cm 2 , respectively. The results indicate 
that the dissociation cross section of cyclohexane is about one or- 
der of magnitude higher when swift heavy cosmic ray analogs were 
employed in comparison with the experiment employing medium- 
mass cosmic rays. 
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(ii) The half-life of cyclohexane (fi/2) in ISM as a result of cos- 
mic ray bombardment for pure c-C 6 H 12 ice is found to be about 
5 x 10 6 years (considering only medium-mass cosmic rays); for 
mixed ice it is about 1 x 10 6 years (considering only heavy-ion cos- 
mic rays). 

(iii) For the mixed ice (H 2 0:NH 3 :c-C 6 H 12 ), after 20 xlO 6 years 
in ISM, almost 20% of the initial cyclohexane was converted into 
CO by heavy cosmic rays, 3% was transformed into OCN~, and 1% 
into C0 2 . Such result suggests that highly hydrogenated hydrocar- 
bons in water-rich grain mantles inside interstellar clouds (or other 
astrophysical ices exposed to medium-mass and heavy cosmic rays) 
can be mostly converted into CO. 

(iv) The results of the radiolysis of simple alkanes such as cy- 
clohexane, in astrophysical ice analogs are consistent with the 
production of unsaturated molecules containing double bond(s), 
triple bond(s) (e.g. 2334 cm -1 ), and/or ring(s) (e.g. 3085 cm 4 ). 
Extrapolating to an astrophysical scenario, the maximum produc- 
tion of alkenes is obtained after about 3-5 xlO 6 years indepen- 
dent of the projectile type and the presence of polar species in 
the ice (e.g. water and ammonia). The formation cross section of 
alkenes, considering the C-H out-of-plane bending mode at 719 
cm 4 (A ~ 8 x 10~ 18 cm molec -1 ), in the pure c-C 6 Hi 2 ice and 
mixed H20:NH3:c-CsHi2 ice were about 1 x 10~ 15 and 1 x 10~ 14 
cm 2 . The maximum value was roughly the same for both experi- 
ments, ~ 10~ 2 alkenes per cyclohexane molecule irradiated by cos- 
mic ray analog. Another consequence of this induced unsaturation 
is the production of molecular hydrogen in the ice, which can be 
released to the gas phase (desorbed) depending on its energy or ice 
temperature. 

(v) The 3.4 fim band of PN IRAS 05341+0852 has large sim- 
ilarity with the band observed for the organic residue (~ 130 K) 
obtained after the ion bombardment of H 2 0:NH 3 :c-C6Hi2 ice. This 
suggests two different issues: first, it puts a constrain in the tem- 
perature of dust grains in the region, and secondly, it also indicates 
that a negligible amount of CO was present in the grains initially. 

(vi) Considering the sharp 3.37 /im absorption peak (C-H asym- 
metric stretching mode in -CH3) we suggest that grains in PPN 
CRL 618, DISM GC IRS 6E, and Murchison carbonaceous me- 
teorite were probably exposed to heavy cosmic rays with fluence 
higher than 3xl0 13 ions cm~ 2 or to some kind of ionizing agents 
that delivered an energy dose higher than ~ 2 x 10 8 J kg' 1 ss 10 15 
eV ng _1 . 

(vii) A comparison between IR spectra from the laboratory and 
YSOs suggests that the initial composition of grains in young stel- 
lar disks may include a mixture of H 2 0, NH 3 , CO (or C0 2 ), alka- 
nes, and CH3OH. 

Finally, the results indicated that cosmic ray bombardment of 
saturated alkanes can be an alternative scenario for the production 
of unsaturated hydrocarbons and possibly aromatic rings (via dehy- 
drogenation processes) in interstellar and protostellar ices. More- 
over, a comparison between the laboratory spectra and infrared ob- 
servations of protostellar ices indicates that saturated hydrocarbons 
such as cyclohexane may be a part of the initial chemical inventory 
of interstellar ices. 
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APPENDIX A: NEW IR FEATURES AFTER RADIOLYSIS 

Table |AT1 lists the new IR features formed after ion bombardment 
of the cyclohexane containing ices at 13 K. 

This paper has been typeset from a TgX/ LTgX file prepared by the 
author. 
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Table Al. Possible assignment of infrared absorption features produced by the radiolysis of the cyclohexane containing ices at 13 K. 



Wavenumber (cm 1 ) 


Attribution and comments 


References 


c-C 6 H, 2 


H 2 0:NH 3 :c-C 6 H 12 






- 


3091 


— v I 1 \ct\J J oil. UC11Z.C11C, Ly V^lUllCACllC, dllvCllCa 


M 61 
U>oj 


3021 


3021" 


— f~"_l-T ^cn^A ctr hpn7Pnp' arrwvt ati r» *? 

— v I 1 \ct\J J BIT. UC11Z.C11C, ill Ul 1 lill 1L . 




- 


3008 


— f 1 _ 1— I / cn^ \ cf i- ol l^f^nt^C 

— v. n ^blj ) bit. cUKtiitB 


n 1 
1 1 1 


- 


2971 


? 




2341 


2337 


C02 


T21 

L^J 


2334 


- 


_ f ^ ~ f ^ _ cvm c.tr of nnn cvm nlKvnfc, 
v — v By in. b u \ji uuu > y 111. ciiiv y ii^b 


t 1 ] 


- 


2258 


CH2-C=N? (C=N str.) 


T71 

L ' J 


- 


2225 


-C=C-H (C=C str of internal alkvnes at 2220 cm -1 )' 

V V 11 - — *. Oil- Ul lllLV^LllUl tllJYV 11^-l} ill ^L.^\J l^lll ) . 


1 1 1 


- 


2166 


OCN- 


pi 


2133 
- 
- 


2137 
2081 
1951 


CO; -C=C-H (C=C str. of terminal alkynes at 2120 cm" 1 )? 
HCN 


rn 

L 1 J 


- 


1713 


7 




1688 


- 


broad C— C str aromatic 1 ^ ring 1 ^ 


Ml 

L 1 J 


1642 


1641 


rirritirl f" 1 — C ctr tirnmntir"'? i"i n cr '? 
Ul UclLl. V — v , all, cllUllltlLIU . , llllg . 


r 1 1 


1584 


1584 


broad NH<:CHiCOCr'? 


141 
t^J 




1436 


1> IT J V, TT7 \J\J ; 


T41 


1392 




L-UULyi, UClliCllC : 


n 1 


1376 


1375 


-CH3 bending sym. t-butyl or dimethyl 


[1] 


1300 


1302 


CH4 (C-H V4 deformation mode) 


[5] 


1137 


1136" 


7 




993 


994 


^C=C^ (C-H bending out plane); cyclohexene; aromatic? 


rn 


967 




? 






954 


? 




948 




? 




916 


916 


f,C=C" (C-H bending out plane); terminal vinyl at 910 cm 4 ? 


[i] 


877 


877° 


7 






845 


f C=C* (C-H bending out plane) 


[i] 


818 




? weak 




810 




? weak 




751 


754 


? broad and strong; aromatic C-H bending out plane; 

CH 2 rocking aliphatic alkanes (present in molecules with less then 4 carbon 
and very strong in molecules with less then 3 carbon) 


[i] 


720 


719" 


^C=C^ bending out of plane aliphatic alkenes; cyclohexene 


[i] 


680 


683 


benzene? (C-H bending in plane) 


[1,6] 


661 




C0 2 bending? 


[2] 


643 


644 


^C=C^ (C-H bending out of plane); cyclohexene 


[1] 



[1] Pavia et al. 2009; [2] Gerakines et al. 1995; [3] Demyk et al. (1998); [4] Holtom et al. 2005; 
[5] Boogert et al. 1997; [6] Ruiterkamp et al. 2005; [7] Larkin 201 1. 
"absent for fluences >3 X 10~ 13 ions cm -2 . 



